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球墨 铸铁 凝固 显 微 组 织 的 元 胞 自动 机 模拟 


Kod nd KAF 
(东南 大 学 江苏 省 先进 金属 材料 高 技术 研究 重点 实验 室 , 南京 211189) 


摘 要 改进 了 前 期 工作 建立 的 多 相 元 胞 自动 机 (multi-phase cellular automaton, MCA) 模 型 , 模拟 以 离异 共 晶 方式 凝固 

墨 铸铁 的 显 微 组 织 演 化 . 在 模型 中 采用 局 部 溶质 平衡 法 计算 石墨 和 奥 氏 体 的 生长 动力 学 , 并 在 石墨 的 生长 模型 中 考虑 石墨 
与 Fe 的 密度 比 . 该 模型 可 以 模拟 出 与 实验 观测 相符 合 的 显 微 组 织 形 貌 . 应 用 该 模型 模拟 分 析 了 石墨 与 奥 氏 体 的 相互 作用 和 
竞争 生长 机 制 , 讨论 了 冷却 速率 对 凝固 结束 时 石墨 球 大 小 和 尺寸 分 布 的 影响 , 将 模拟 结果 与 实验 结果 进行 了 比较 . 结果 
表明 : 奥 氏 体 的 析出 促进 邻近 石墨 在 液 相 中 的 生长 ; 奥 氏 体 和 石墨 两 相 的 生长 受 C 扩 散 控制 ; 当 石 墨 被 奥 氏 体 包 围 后 ， 
生长 速度 减 慢 . 此 外 , 随 着 冷却 速率 的 增 大 , 凝固 时 间 缩 短 , 石墨 球 平均 半径 减 小 , 不 同 冷 速 条 件 下 石墨 球 尺寸 分 布 的 变 
化 规律 与 实验 结果 吻合 较 好 . 
关键 词 球墨 铸铁 , 凝固 , 离异 共 晶 , 元 胞 自动 机 , 显 微 组 织 模拟 
中 图 法 分 类 号 TG161 文献 标识 码 A 


ER 


(er 


文章 编号 0412-1961(2015)02-0148-11 


SIMULATION OF SOLIDIFICATION MICROSTRUC- 
TURE OF SPHEROIDAL GRAPHITE CAST IRON 
USING A CELLULAR AUTOMATON METHOD 


ZHANG Lei, ZHAO Honglei, ZHU Mingfang 
Jiangsu Key Laboratory for Advanced Metallic Materials, Southeast University, Nanjing 211189 
Correspondent: ZHU Mingfang, professor, Tel: (025)83793355, E-mail: zhumf(üseu.edu.cn 


Supported by National Natural Science Foundation of China (No.51371051) 
Manuscript received 2014-06-17, in revised form 2014-12-10 


ABSTRACT Spheroidal graphite (SG) cast iron is characterized by the presence of spherical graphite nodules 
distributed in the metallic matrix. The performance of castings is primarily dependent on the solidification micro- 
structures. In this work, a two dimensional (2D) multi-phase cellular automaton (MCA) model previously proposed 
by the present authors is improved to simulate the microstructure evolution of SG cast iron during divorced eutec- 
tic solidification. The present model adopts a local solutal equilibrium approach to calculate the driving force for 
the growth of both graphite and austenite phases. The density difference between iron and graphite is also taken in- 
to account. The diffusion of solute in the simulation domain is calculated using a finite difference method (FDM). 
The present model is applied to simulate the evolution of microstructure and carbon concentration field during so- 
lidification for hypereutectic SG cast irons. The results show that the present model can reasonably describe the typ- 
ical features of divorced eutectic solidification, involving the independent nucleation and growth of primary graph- 
ite and austenite dendrites in liquid, the competitive growth of adjacent graphite nodules, engulfment of graphite 
nodules by austenite dendrites, the isotropic growth of the austenite shells that envelop the graphite nodules, the 
austenite to graphite eutectic phase transition controlled by carbon diffusion through the solid austenite shell, and 
multiple graphite nodules encapsulated in each austenite grain at the end of eutectic solidification. The simulated 


volume fraction and average diameter for graphite nodules are compared reasonably well with the experimental da- 
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ta and level rule calculation. The interactive and competitive growth behavior between austenite dendrites and 


graphite nodules is studied in detail. It is found that the growth of a graphite nodule is promoted by the approach- 


ing austenite. However, after embedded by an austenite dendrite, the growth velocity of graphite decreases rapidly 


because of lower carbon diffusivity in austenite than that in liquid. In addition, the effect of cooling rate on the size 


of graphite nodules is also investigated. The results show that with cooling rate increasing, the size distribution of 


graphite nodules varies from two peaks to one peak, and the average diameter of nodules decreases. The simulation 


results compare reasonably well with the experimental data reported in literature, demonstrating the validity of the 


present model. 


KEY WORDS spheroidal graphite cast iron, solidification, divorced eutectic, cellular automaton, microstructure 
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球墨 铸铁 中 的 石墨 为 球形 , 与 片 层 状 尖 锐 石 黑 
相 比 , 不 仅 应 力 集 中 小 , 而 且 在 同样 的 石墨 体积 下 ， 
球状 石墨 的 表面 积 最 小 , 因而 对 基体 的 割裂 作用 也 
降 到 最 小 . 球墨 铸铁 自问 世 以 来 , 由 于 在 力学 性 能 、 
铸造 性 能 和 价格 方面 的 优势 , 在 机 械 、 冶 金 、 汽 车 、 
矿山 等 工业 部 门 得 到 了 广泛 的 应 用 ". 

球墨 铸铁 的 宏观 力学 性 能 与 其 在 凝固 过 程 中 
所 形成 显 微 组 织 密切 相关 , 包括 石墨 球形 状 、 大 小 、 


根据 平衡 相 图 确定 界面 浓度 , 没有 耦合 实时 浓度 场 
计算 . 
Burbelco 等 (应 用 连续 形 核 模 型 , 将 动力 学 过 
冷 作为 生长 驱动 力 , 模拟 了 球墨 铸铁 凝固 过 程 奥 氏 
体 和 石墨 的 形 核 与 生长 . 模拟 结果 显示 了 与 实验 观 
测 相 符 的 不 规则 圆 形 石 墨 形 貌 , 以 及 凝固 过 程 的 冷 
却 曲 线 . 但 该 模型 没有 针对 冷却 速率 等 工艺 因素 对 
凝固 显 微 组 织 的 影响 规律 开展 研究 , 并 且 模 型 中 没 
考虑 对 于 凝固 组 织 的 形成 有 着 重要 影响 的 石墨 


c— 


Hu 


尺寸 分 布 以 及 基体 组 织 的 形 貌 . 工业 上 应 用 较 广 泛 
的 是 过 共 唱 球墨 铸铁 . 在 其 凝固 过 程 中 , 当 温 度 下 
降 到 石墨 的 液 相 线 温 度 时 , 初生 石墨 相 在 液 相 中 形 
核 并 生长 , 直到 共 唱 温度 , 奥 氏 体 开始 形 核 . 奥 氏 体 
晶 核 与 石墨 接触 之 前 在 液 相 中 以 枝 晶 形式 生长 中 . 
石墨 球 与 奥 氏 体 发 生 接 触 之 后 , 石墨 球 迅速 被 奥 氏 
体 壳 包围 , 并 在 奥 氏 体 中 继续 生长 . 由 于 C 在 奥 氏 体 
中 扩散 较 慢 , 石墨 被 奥 氏 体 包 围 后 生长 缓慢 . 同时 ， 
奥 氏 体 壳 在 液 相 中 继续 生长 , 直到 奥 氏 体 壳 相互 接 
fih, 最 终 完 成 凝固 过 程 . 

目前 , 计算 机 模拟 方法 已 经 发 展 成 为 研究 凝固 
过 程 中 显 微 组 织 演变 机 理 的 一 种 有 效 的 方法 , 其 中 
元 胞 自动 机 (cellular automaton, CA) 方 法 已 成 为 该 
领域 研究 的 热点 . 国内 外 学 者 们 提出 了 各 种 二 维和 
三 维 的 CA 模型 , 对 纯 扩散 和 对 流 作 用 下 的 枝 晶 生 
长 、 微 观 偏 析 、 显 微 气孔 的 生长 等 凝固 组 织 进 行 模 
WR, Charbon 和 Rappaz09 最 早 将 CA 方法 应 用 
到 球墨 铸铁 凝固 过 程 , 可 视 化 地 研究 凝固 过 程 中 显 
微 组 织 的 演变 , 他 们 基于 C 扩 散 控制 生长 , 应 用 经 典 
模型 模拟 球 铁 凝固 过 程 . 其 模拟 结果 显示 , 每 个 共 
品 唱 粒 中 只 包含 一 个 石墨 球 , 这 与 实验 观测 结果 不 
符 . Ruxanda 等 "在 此 模型 基础 上 考虑 了 初生 奥 氏 
体 的 影响 , 并 假设 石墨 与 奥 氏 体 接触 之 后 才 开 始 生 
K, 对 亚 共 晶 球 墨 铸铁 的 凝固 过 程 进行 模拟 , 模拟 
结果 中 每 个 奥 氏 体 唱 粒 包 含 多 个 石墨 . 但 该 模型 没 
有 考虑 石墨 在 液 相 中 的 生长 , 且 计 算 生 长 动力 学 时 


与 Fe 两 相 的 密度 比 . 

本 课题 组 前 期 工作 ”依据 Zhu 和 Stefanescu"” 
提出 的 基于 局 部 溶质 平衡 的 模型 , 以 C 扩 散 作 为 奥 
氏 体 和 石墨 生长 的 驱动 力 , 提出 了 针对 球墨 铸铁 凝 
固 过 程 中 球状 石墨 和 奥 氏 体 枝 唱 的 形 核 和 生长 的 
多 相 CA 模型 . 由 于 在 石墨 生长 的 模型 中 人 为 地 将 
石墨 界面 的 C 成 分 进行 平均 化 处 理 , 模拟 结果 中 石 
墨 呈 现 规整 的 球形 . 模型 中 没有 考虑 到 石墨 与 Fe 的 
密度 比 , 使 得 模拟 结果 中 石墨 体积 分 数 和 平均 半径 
均 低 于 实验 测量 的 数据 . 

本 工作 在 前 期 工作 的 基础 上 , 将 多 相 元 胞 自动 
机 (multi-phase cellular automaton, MCA) 模 型 作 了 改 
进 , 在 石墨 的 生长 动力 学 中 采用 石墨 界面 实际 的 局 
部 C 浓 度 进行 计算 , 并 考虑 了 石墨 与 Fe 的 密度 比 ， 
模拟 了 过 共 蝇 球墨 铸铁 凝固 过 程 中 显 微 组 织 形 貌 
与 区 域 浓度 场 的 变化 , 分 析 了 石墨 与 奥 氏 体 的 相互 
作用 和 竞争 生长 机 理 . 应 用 该 模型 模拟 研究 了 冷却 
速率 对 于 凝固 显 微 组 织 中 石墨 球 的 尺寸 分 布 以 及 
平均 半径 的 影响 规律 . 

1 模型 与 算法 
1.1 模型 描述 

本 工作 主要 研究 过 共 唱 球墨 铸铁 以 离异 共 唱 
方式 凝固 时 球状 石墨 与 奥 氏 体 两 相 之 间 的 相互 作 
与 竞争 机 制 以 及 浓度 场 的 变化 规律 . 将 二 维 的 计 
算 区 域 剖 分 成 均匀 的 正方 形 网 格 , 在 液 相 、 奥 氏 体 
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相 和 石墨 相 三 相 系 统 中 , 每 个 网 格 被 赋予 浓度 、 温 。 长 界面 前 沿 存在 浓度 梯度 , 引起 溶质 C 的 扩散 . 因 
度 . 3 个 相 ( 液 相 、 奥 氏 体 和 石墨 ) 的 体积 分 数 、 奥 氏 体 e, 考虑 了 溶质 再 分 配 的 溶质 传输 控制 方程 为 : 


枝 晶 的 择优 取向 、 石 墨 球 标号 等 变量 , 并 标明 元 胞 
的 状态 : 液 相 (=1), 8 SR TH (1), K REHE), 
石墨/ 液 相 界 面 (0< 包 二 1, 50), 石墨 / 奥 氏 体 界 面 (0< 
fa 1, f=0), 奥 氏 体 / 液 相 界面 (0<f; 达 1, fa=0), 石墨/ 
奥 氏 体 / 液 相 界 面 (6>0, f>0, fofi 1), 其 中 fi ,i 和 
JJ 分别 代 表 液 相 、 石 墨 和 奥 氏 体 的 相 分 数 . 
模拟 开始 时 , 首先 将 区 域内 网 格 初始 化 , 使 所 
有 网 格 具 有 相同 的 成 分 和 温度 , 本 工作 中 温度 计算 
从 初始 成 分 所 对 应 的 液 相 线 温 度 开始 . SET EE h 
球墨 铸铁 , 当 区 域 温度 低 于 石墨 液 相 线 温度 时 , 石 
墨 从 液 相 中 析出 , 并 以 近似 球状 长 大 . 当 温 度 降 至 
dt dia EE, 奥 氏 体 析 出 并 和 石墨 相 一 起 以 离异 共 唱 
方式 生长 . 石墨 和 奥 氏 体 的 形 核 密度 根据 文献 [18] 
的 形 核 模 型 或 相关 的 实验 数据 确定 . 石墨 和 奥 氏 体 
形 核 后 位 置 固定 , 不 考虑 流 场 而 造成 石墨 的 上 浮 作 
用 . 在 凝固 过 程 中 , 当 石墨 与 奥 氏 体 相 接触 时 , 石墨 
迅速 被 奥 氏 体 壳 包围 . 随后 , 石墨 通过 周围 奥 氏 体 
过 吸收 溶质 C 继 续 生 长 , 直到 共 析 温度 模拟 结束 . 奥 
氏 体 生长 过 程 中 不 断 向 周围 液 相 中 排出 溶质 C, 直 
至 奥 氏 体 壳 相互 接触 , 完成 凝固 . 本 模型 不 考虑 流 
动 对 溶质 传输 的 影响 , C 溶 质 的 传输 通过 纯 扩 散 实 
现 . 假设 计算 区 域内 温度 均匀 分 布 , 并 以 一 定 的 冷 
却 速率 冷却 . 凝固 过 程 中 的 冷却 速率 受到 散热 速 
率 、 石 墨 与 奥 氏 体 生长 过 程 中 释放 的 潜 热 两 方面 作 
用 的 共同 影响 . 
本 工作 采用 基于 局 部 溶质 平衡 模型 "计算 石墨 
与 奥 氏 体 的 生长 , 两 相 的 生长 驱动 力 为 界面 平衡 成 
分 与 实际 成 分 的 差 值 相关 的 函数 . 其 中 局 部 界面 平 
衡 成 分 根据 界面 局 部 温度 和 界面 曲率 计算 得 到 , 而 
实际 成 分 由 浓度 场 计 算 获得 . 此 外 , 本 模型 在 计算 
两 相 的 生长 时 考虑 了 石墨 与 Fe 的 密度 比 . 
1.2 浓度 场 和 温度 场 计算 
假设 在 奥 氏 体 / 液 相 ` 石 墨 / 液 相 和 石墨 / 奥 氏 体 
界面 上 的 溶质 平衡 . 一 方面 , 奥 氏 体 在 生长 过 程 中 
向 周围 液 相 排 出 溶质 , 固 / 液 界面 的 溶质 再 分 配 的 计 
算 如 下 : 


| 回 


GC =k, C, (1) 
式 中 ,无 为 奥 氏 体 / 液 相 界面 溶质 平衡 分 配 系数 ，C， 
和 C, 分 别 为 奥 氏 体 / 液 相 界面 处 奥 氏 体 和 液 相 的 成 
分 . 男 一 方面 , 石墨 在 生长 过 程 中 从 周围 的 液 相 或 
者 奥 氏 体 相 中 吸收 溶质 C. 两 者 的 共同 作用 使 得 生 


=Y(DYCJ+CIU- DEA -C, jm (2) 
式 中 ,1! 是 时 间 , C; 是 成 分 , D; 是 扩散 系数 , 下 标 i 表 示 
奥 氏 体 相 y 或 者 液 相 L; C, 为 石墨 生长 时 石墨 / 奥 氏 
体 界面 中 奥 氏 体 相 的 成 分 , 或 者 石墨 / 液 相 界 面 和 石 
墨 / 奥 氏 体 / 液 相 界面 中 的 液 相 成 分 ; C6 为 石墨 的 C 
含量 (Ce=100%). 等 式 右边 第 二 项 表示 奥 氏 体 相 生 
长 时 在 奥 氏 体 / 液 相 或 奥 氏 体 /石墨 / 液 相 界 面 产生 的 
溶质 再 分 配 . 其 处 理 方 法 为 : 将 奥 氏 体 生 长 所 释放 
的 溶质 添加 到 相同 界面 网 格 的 剩余 液 相 中 . 如 果 该 
界面 网 格 在 该 时 间 步 时 已 完全 凝固 , 则 将 所 释放 的 
溶质 添加 到 与 其 相 邻 的 液 相 网 格 中 . 等 式 右边 第 三 
项 为 石墨 生长 时 在 石墨 / 液 相 ` 石墨/ 奥 氏 体 或 石墨 / 
奥 氏 体 / 液 相 界 面 产生 的 溶质 再 分 配 . 其 处 理 方 法 
23: 将 石墨 生长 所 吸收 的 溶质 从 相同 界面 网 格 的 剩 
余 液 相 或 奥 氏 体 相 中 减 去 相应 的 溶质 量 . 如 果 该 界 
面 网 格 在 该 时 间 步 时 已 完全 转变 成 石墨 , 则 将 所 吸 
收 的 溶质 从 与 其 相 邻 的 液 相 或 奥 氏 体 网 格 中 减 去 
相应 的 溶质 量 

在 模拟 体系 中 包含 了 液 相 、 石墨 和 奥 氏 体 3 个 
不 同 的 相 , 以 及 奥 氏 体 / 液 相 、` 石 墨 / 液 相 、 石 墨 / 奥 氏 
体 和 石墨 / 奥 氏 体 / 液 相 4 种 不 同 的 界面 . 在 不 同 的 相 
和 界面 处 的 C 成 分 和 C 扩 散 系数 差别 是 很 大 的 . 为 
了 有 效 合理 地 处 理 C 成 分 和 扩散 系数 在 界面 处 的 不 
连续 问题 , 在 处 理 溶质 扩散 时 , 采用 等 效 成 分 C. 和 
等 效 扩 散 系 数 D. 求 解 式 (2). 例如 , 液 相 和 奥 氏 体 的 

等 效 成 分 分 别 为 C=Ct 和 C=C/k; 液 相 和 奥 氏 体 的 

等 效 扩 散 系 数 分 别 为 D=Dr 和 D=D; 石墨 相 中 及 
其 与 之 相 邻 各 相 之 间 的 等 效 扩散 系数 为 C=0. 这 
是 因为 石墨 相 的 平衡 成 分 是 100%C. 虽然 与 其 周 
围 的 液 相 或 奥 氏 体 相 之 间 的 浓度 相差 很 大 , 但 根 
据 相 平衡 原理 , 它们 之 间 的 化 学 位 是 相等 的 . 因 
此 , 在 与 石墨 相 邻 的 界面 处 不 会 发 生 C 从 石墨 到 
液 相 或 奥 氏 体 的 扩散 . 计算 区 域 中 网 格 的 种 类 和 
不 同 种 类 网 格 中 的 平衡 浓度 与 扩散 系数 的 详细 确 
定 方法 如 文献 [19] 所 述 . 

采用 显示 格式 的 有 限 差 分 法 和 无 扩散 边界 条 
件 对 式 (2) 进 行 离散 和 求解 , 稳 态 的 时 间 步 长 At 由 C 
在 液 相 中 的 扩散 系数 和 网 格 尺 寸 Ax 确定 : 

At = Ax /(4.5D,) 3) 

在 凝固 过 程 中 , 一 方面 铸件 通过 铸 型 与 外 界 发 

生 热 交换 , 引起 体系 温度 的 下 降 , 为 凝固 提供 驱动 
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73; 男 一 方面 , 凝固 过 程 伴 随 着 相 变 潜 热 的 释放 , 引 


起 温度 的 回升 . 这 两 方面 的 共同 作用 使 得 体系 的 温 
度 随 着 凝固 时 间 而 变化 . 此 外 , 不 同 温度 下 和 铸件 内 
各 点 的 位 置 和 状态 也 会 存在 差异 . 不 过 , 由 于 本 工 
作 的 模拟 区 域 比较 小 , 并 且 Fe-C 合 金 的 热 扩 散 系 数 
比 溶质 C 在 液 相 中 的 扩散 系数 大 约 3 个 数量 级 , 所 
以 , 本 模型 假设 在 一 个 微观 区 域内 热 扩 散 能 够 充分 
进行 , 计算 区 域内 温度 均匀 , 区 域 中 温度 受 散 热 速 
率 与 潜 热 释放 两 方面 共同 作用 , 区 域 中 温度 了 的 计 
算 公式 为 : 


oT $4 Ly Uf, us Of (4) 
ðt Pre€, Ob Pre, Ôt 


式 中 , 了 表示 散热 速率 ， 指 的 是 铸件 通过 铸 开 与 外 
界 发 生 热 交 换 , 引起 体系 温度 下 降 的 速率 , 可 根据 
发 生 相 变 前 和 相 变 结束 后 的 冷却 曲线 拟 合 确 定 ; pi 
为 Fe 的 密度 ; cs 为 Fe 的 比热容 ; 六 ec M Lpy DAA 
石墨 与 Fe 的 体积 潜 热 . 
1.3 石墨 生长 模型 
球墨 铸铁 凝固 过 程 中 , 石墨 的 生长 分 为 在 液 相 
中 生长 和 被 奥 氏 体 包围 后 在 奥 氏 体 壳 中 生长 2 个 阶 
段 . 如 前 文 所 述 , 本 工作 采用 基于 局 部 溶质 平衡 模 
型 中 计算 石墨 的 生长 , 生长 驱动 力 为 局 部 界面 平衡 
成 分 与 实际 成 分 之 差 . 当 石 墨 在 液 相 中 生长 时 , 需 
要 考虑 界面 曲率 的 影响 , 石墨 / 液 相 界面 的 平衡 成 分 


Con 按 下 式 进 行 计算 : 


q 
(i E C E: T'- Ex xL, eem 
Gr/L 0 


式 中 , 7’ 为 界面 温度 ,76% 为 初始 成 分 Co 所 对 应 的 
液 相 界 面 的 界面 曲率 , To 7348 581] Gibbs- Thomson 
系数 . 在 本 工作 的 2D 模型 中 , 石墨 界面 平均 曲率 由 


下 式 计算 : 
Kon. = 一 志 一 0 
ANA (6) 


RP, 大 为 每 个 元 胞 的 面积 ，Y 太 , 为 属于 同一 个 石 


墨 的 石墨 体积 分 数 总 和 
根据 界面 溶质 守恒 , 并 考虑 石墨 和 Fe 之 间 的 密 
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AF, pu 为 石墨 的 密度 ; Cua 为 石墨 / 液 相 界 面 的 
成 分 , 通过 浓度 场 计算 获得 ; Ge 为 与 邻 位 网 格 状态 
相关 的 几何 因子 , 由 下 式 确定 ; 
Go, -uini (xs. l Lis] 
I gl o XI 
s.s zn Vie 
AP, S 为 最 邻近 的 4 个 网 格 的 状态 参数 ; S" 为 次 近 
邻 的 4 个 网 格 的 状态 参数 . 几何 因子 Ce 考虑 了 邻 位 
的 石墨 相 网 格 数 对 该 网 格 的 石墨 相 增 长 的 影响 , 同 
时 还 考虑 了 邻 位 网 格 的 影响 程度 与 网 格 间距 有 关 ， 
即 最 近邻 网 格 比 次 近邻 的 网 格 的 影响 更 大 . 
当 石墨 在 奥 氏 体 中 生长 时 , 只 需 将 式 (5) 中 TO% 
FI mis. 分 别 换 成 石墨 / 奥 氏 体 界 面 的 平衡 温度 Te 
和 石墨 的 固 相 线 斜率 msov 并 且 忽 略 界面 曲率 的 作 
], 其 余 均 与 石墨 在 液 相 中 生长 的 计算 方法 相同 . 
因此 , 石墨/ 奥 氏 体 界 面 的 平衡 成 分 Cey 为 : 


(8) 


pup 
Goss ——n T (9) 


与 式 (7) 的 计算 方法 类 似 , 一 个 时 间 步 长 内 , 厂 
黑 / 奥 氏 体 界面 的 石墨 体积 分 数 增 量 Afi 为 : 
(Ces Ces ) 
Ag cg POR COM 10 
a RET (10) 
RH, C, 为 石墨/ 奥 氏 体 界面 的 奥 氏 体 的 局 部 平 
均 成 分 . 对 于 石墨 / 奥 氏 体 的 界面 元 胞 其 局 部 平均 
成 分 可 按 下 式 计算 : 
200 GS Pe) 
Cop (I) = E (11) 
RP, CU) 为 界面 元 胞 7 的 奥 氏 体 相 成 分 ， 
Ci,(J) 为 元 胞 7 的 邻 位 界面 元 胞 的 奥 氏 体 相 成 
J=1 
分 之 和 , N 为 邻 位 界面 元 胞 的 个 数 . 局 部 的 奥 氏 体 成 
分 C, U) 和 C U) 由 浓度 场 计算 获得 . 
当 石墨 / 液 相 界 面 或 者 石墨/ 奥 氏 体 界面 元 胞 的 


石墨 体积 分 数 达到 1 时 , 该 元 胞 状态 即 从 石墨 / 液 相 
界面 或 者 石墨 / 奥 氏 体 相 界面 元 胞 转变 成 石墨 , 并 且 


度 差 , 在 一 个 时 间 步 长 内 , 石墨 / 液 相 界 面 网 格 的 石 
墨 体积 分 数 增 量 A 太 .为 : 


E 
Afo = GaP t a 2) (7) 
Pe (Cs Com) 


REJA E29 E SACRI PR TRU rc LEG 89/9. E 
相 界 面 元 胞 , 从 而 使 该 石墨 得 以 继续 生长 . 
1.4 奥 氏 体 生长 模型 

本 工作 同样 采用 局 部 溶质 平衡 法 "计算 奥 氏 体 
枝 晶 的 生长 动力 学 , 考虑 了 枝 晶 的 择优 生长 方向 ， 
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界面 曲率 变化 以 及 界面 能 和 生长 动力 学 的 各 向 异 
性 . 枝 品 生长 的 驱动 力 是 极品 / 液 相 界面 实际 成 分 与 
平衡 成 分 之 差 , 在 液 相 中 生长 时 枝 唱 界面 平衡 成 分 
的 Cy, 计算 方法 如 下 : 

Fem. VES fi(0,0,) 


mL 


CA t 


YAL 


12 
aY m, l ) 


式 中 , T 为 奥 氏 体 / 液 相 界 面 的 平衡 温度 ; mL, 为 奥 
氏 体 的 液 相 线 斜率 ; 九 为 奥 氏 体 的 Gibbs-Thomson % 
数 ; K, 为 奥 氏 体 / 液 相 界 面 曲率 , 通过 计数 格子 法 中 计 
算 获 得 ; f(0.0,) 为 界面 能 各 向 异性 函数 , 由 下 式 进 
行 计 算 : 


fi(0,00)=1-6, cos|4(6 - 6,)] 


" f, /ox (13) 
(af, ax? (af, ay)" 


式 中 , 6. 是 界面 能 各 向 异性 强度 , & 是 枝 唱 的 择优 生 
长 方向 . 

根据 界面 溶质 平衡 , 在 一 个 时 间 步 长 内 , 奥 氏 
体 / 液 相 界 面 元 胞 的 奥 氏 体 体 积分 数 增 量 Aj 为 : 
(655€) 
(c7. 3) 
式 中 , Cj 为 根据 式 (12) 计 算 的 奥 氏 体 / 液 相 界面 的 
平衡 成 分 ; Cn 为 由 浓度 场 计算 获得 的 奥 氏 体 / 液 相 
界面 的 实际 成 分 ; G, 为 几何 因子 , 计算 方法 与 式 (8) 
相似 , 只 需 把 式 (8) 中 的 石墨 体积 分 数 换 成 奥 氏 体 体 
积分 数 ; 万 (6,0) 为 界面 动力 学 各 向 异性 函数 , 计算 
方法 如 下 : 

f(0,00)=1+6, cos[4(0 - 6.) (15) 

AP, 6, 是 界面 动力 学 各 向 异性 强度 . 

当 奥 氏 体 / 液 相 界面 元 胞 的 奥 氏 体 体积 分 数 达 
到 1 时 , 该 元 胞 的 状态 由 奥 氏 体 / 液 相 界 面 元 胞 转变 
成 奥 氏 体 元 胞 , 并 捕捉 周围 的 液 相 元 胞 为 新 的 奥 氏 
体 / 液 相 界面 元 胞 , 使 得 奥 氏 体 在 下 一 个 时 间 步 长 内 
继续 生长 . 当 石墨 在 奥 氏 体 中 生长 时 , 一 个 时 间 步 
KA, 石墨 / 奥 氏 体 界面 元 胞 中 石墨 体积 分 数 增 量 等 
于 奥 氏 体 相 减少 的 量 . 
1.5 石墨 与 奥 氏 体 枝 晶 生长 模型 的 耦合 
将 石墨 与 奥 氏 体 相 的 形 核 、 生 长 与 浓度 场 和 温 
度 场 计算 相 耦 合 的 方法 如 下 : 模拟 开始 时 , 首先 把 
计算 区 域内 所 有 网 格 初 始 化 , 使 其 具有 相同 的 初始 
浓度 和 温度 . 然后 根据 散热 速率 计算 区 域 以 一 定 的 


Af, = G, f. (6,6, (14) 
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冷却 速率 冷却 , 根据 不 同 的 实时 温度 判断 是 否 符合 
石墨 与 奥 氏 体 的 形 核 条 件 , 如 果 符 合 , 石墨 或 奥 氏 
体 开始 形 核 , 赋予 这 些 形 核 元 胞 相应 的 状态 值 . 这 
时 存在 多 种 不 同 的 界面 , 当 奥 氏 体 相 生长 时 , 根据 
式 (12)~(15) 计 算 奥 氏 体 相 分 数 的 增 量 . 石墨 在 液 相 
中 和 奥 氏 体 相 中 的 生长 分 别 根 据 式 (5)~(8) 和 (9)~ 
(11) 计 算 . C 溶 质 扩 散 和 由 于 相 变 造成 的 溶质 再 分 配 
根据 式 (1) 和 和 (2) 进行 计算 . 如 前 所 述 , 本 模型 考虑 了 
石墨 与 Fe 的 密度 比 . 为 了 使 整个 计算 区 域 中 的 溶质 
FE, 在 石墨 生长 后 进行 溶质 再 分 配 时 , 用 等 效 的 
石墨 C 溶 度 ( C5.=100%x (pc. /p,. ) ) 进 行 计算 . 由 于 
CA 方法 中 对 于 溶质 再 分 配 的 处 理会 造成 局 部 溶质 
分 配 的 误差 , 每 个 时 间 步 长 对 区 域内 溶质 的 总 量 进 
行 统计 和 补偿 , 随后 计算 整个 区 域内 溶质 C 的 扩散 . 
如 前 文 所 述 , 散热 速率 的 作用 使 区 域 温 度 下 降 , 而 
石墨 和 奥 氏 体 的 析出 伴随 着 潜 热 的 释放 , 引起 区 域 
内 温度 升 高 , 两 者 的 共同 作用 对 区 域 温度 的 影响 用 
式 (4) 进 行 计算 . 温度 场 和 浓度 场 的 变化 又 会 影响 石 
墨 与 奥 氏 体 的 形 核 和 生长 , 并 开始 新 一 轮 循 环 . 所 
以 , 本 模型 包含 了 一 个 耦合 浓度 场 和 温度 场 计算 、 
石墨 与 奥 开 体形 核 、 生 长 以 及 相互 作用 过 程 的 球墨 
铸铁 离异 共 吕 凝固 显 微 组 织 演变 的 物理 机 制 . 本 工 
作 计 算 所 采用 的 Fe-C 合 金 的 物性 参数 见 表 109. 
2 结果 与 讨论 
2.1 过 共 晶 球墨 铸铁 凝固 显 微 组 织 演变 

在 本 工作 中 , 为 了 将 部 分 模拟 结果 和 实验 结果 
相 比 较 , FH C 当量 来 考虑 Si 元 素 的 作用 . C 当量 用 式 
C"=%C+0.3%Si "进行 计算 . 图 1 为 模拟 的 C 当量 
在 C=4.71% (质量 分 数 ) 的 球墨 铸铁 凝固 显 微 组 织 
和 浓度 场 的 演化 过 程 . 计算 区 域 由 200x200 个 网 格 
组 成 , 网 格 尺寸 为 2 pm. 为 了 与 实验 照片 进行 更 好 
地 比较 , 针对 图 1 的 模拟 时 , 区 域内 部 分 石墨 与 奥 
氏 体 的 形 核 位 置 人 为 设 定 . 初始 温度 为 G, 成 分 所 
对 应 的 石墨 液 相 线 温度 1357 "C. 当 Fe 液 温度 低 于 
石墨 的 液 相 线 温度 后 , 随 着 温度 的 下 降 石墨 不 断 从 
液 相 中 析出 , 当 温 度 达 到 共 兄 温度 时 , 奥 氏 体 开 
始 析 出 . Fe 液 在 温度 高 于 (ZL-10) C 时 的 散热 速率 
X —4.6exp(- 0.00312), 在 该 温度 以 下 的 散热 速率 为 
1 'C/s. 温度 场 受 到 散热 速率 和 石墨 与 奥 氏 体 凝 固 
过 程 中 放出 的 潜 热 共同 影响 . 图 1 中 的 数字 为 局 部 
C 浓 度 . 可 以 看 出 , 凝固 开始 时 , 石墨 作为 初生 相 在 
液 相 中 形 核 生 长 . 石墨 的 生长 吸收 溶质 C. 所 以 , 使 
区 域内 C 浓 度 低 于 初始 成 分 C=4.71%, 并 且 石 墨 周 


I 


a a 
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表 1 本 工作 计算 所 采用 的 物性 参数 ” 


Table 1 Physical properties used in the simulation 


[19] 


Parameter Value Unit 
Partition coefficient k, 0.49 
Liquidus slope of graphite mior 470 K- 95" 
Solidus slope of graphite mso 284.8 K: 95" 
Liquidus slope of austenite m, 一 90 K: 947 
Eutectic temperature T. 1153.744.865x94Si C 
Eutectic composition C. 4.28 % 
Gibbs-Thomson coefficient of graphite To: 1.9x107 m: K 
Gibbs-Thomson coefficient of austenite 7, 3.7x107 mK 
Specific heat of liquid c, 880 本 kg … K` 
Volume latent heat of austenite Ly, 1.82x10* J m^ 
Volume latent heat of graphite Lyc. 3.60x10* Jm? 
1.5x10" — (T»(14.—200) 'C) 
Diffusion coefficient in austenite D, 1.5x10"exp(-142100/8.3147) m^ s" 
(T&(14-200) 'C) 
Diffusion coefficient in liquid D; 5.0x10^? m° s” 
Degree of the surface energy anisotropy à 0.6 
Degree of the kinetic anisotropy & 0.6 
Density of graphite o. 2100 kg: m? 
Density of iron pre 7930 kg: m? 


1 不 同 温 度 下 初始 成 分 C=4.71% 时 过 3 
Fig.1 Simulated and experimental morphologies for a hypereutectic s 
centration Co=4.71% at 7=1253 °C with f=1% (a), T-1163 


Ho ERAS PEE IE RZE 
pheroidal graphite (SG) cast iron with an initial con- 
C with f=3% (b), T-1165 "C with f=6% (c), T= 


1156 'C with f=35% (d), T=1117 'C with f=91% (e), T-1123 "C with f=97% (f), T2740 °C with f=100% (g) 


and experiment result of Fe-3.8696C-2.8296 Si" (h) (Numbers in the figures show the local carbon concentration, 


and f is the total solid fraction) 


LN 
IP 
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围 C 浓 度 略 低 于 远离 石墨 的 区 域 (图 la). 当 温度 降 
至 共 唱 温度 时 , 奥 氏 体 开始 形 核 , 此 时 区 域内 的 浓 
度 接近 共 唱 成 分 Cur=4.28%, 并 且 奥 氏 体 附 近 的 液 
相 中 C 成 分 略 高 , 而 石墨 附近 液 相 中 的 C 成 分 略 低 
(图 1b). 当 石 墨 与 奥 氏 体 接触 时 , 石墨 迅速 被 奥 氏 体 
所 包围 (图 1c), 之 后 奥 氏 体 逐 渐 从 枝 晶 状 转变 为 团 
REK, 溶质 C 从 液 相 经 过 奥 氏 体 过 传输 到 石墨 周 
围 , 使 石墨 在 奥 氏 体内 继续 生长 . 当 大 多 数 石墨 被 
奥 氏 体 包裹 后 , 液 相 中 的 C 浓 度 逐 渐 升 高 (图 1d~e). 
这 是 因为 奥 氏 体 生 长 过 程 中 排出 溶质 C, 而 C 在 奥 
氏 体 内 的 扩散 系数 远 小 于 其 在 液 相 中 的 扩散 系数 . 
所 以 , 奥 氏 体 排出 的 溶质 大 多 进入 液 相 中 , 使 得 液 
相 浓度 升 高 . 凝固 末期 ， 奥 开 体 这 相互 接触 , 交 [i] £i 
L 每 个 奥 氏 体 晶 粒 中 包含 多 个 石墨 . 之 后 温度 继 

续 下 降 , 由 于 奥 氏 体内 浓度 不 均匀 , 石墨 通过 吸收 
IS PA C A EKK, 并 且 石 墨 / 奥 氏 体 界 面 的 
浓度 低 于 远离 石墨 区 域内 浓度 (图 1g). 对 比 740 "C 
时 的 模拟 结果 (图 1g) 与 实验 照片 (图 1b) 呈 可 以 看 
出 , 本 模型 可 以 模拟 出 与 实验 观察 比较 吻合 的 球 铁 
显 微 组 织 形 貌 . 
图 2 是 凝固 过 程 中 国 相 分 数 、 奥 氏 体 和 石墨 体 
积分 数 随 温度 变化 的 曲线 . 可 以 看 到 , 共 晶 温度 之 
前 , 石墨 在 液 相 中 生长 , 石墨 体积 分 数 增长 . 当 温 度 
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图 2 固 相 、 奥 氏 体 和 石墨 体积 分 数 随 温度 的 变化 


Fig.2 Volume fractions of solid, austenite and graphite 


vary with temperature for a hypereutectic SG cast 
iron with C= 4.7196 


ChinaXiv 合 作 期 刊 


学 报 第 51 郑 


达到 共 唱 温度 时 , 奥 氏 体形 成 , 固 相 分 数 、 奥 氏 体 和 
石墨 体积 分 数 迅 速 增长 , 说 明 凝 固 过 程 主要 在 共 唱 
阶段 完成 . 当 固 相 体积 分 数 达 到 100% 后 , 石墨 体积 
分 数 继续 缓慢 增加 , 奥 氏 体 体 积分 数 减少 , 说 明 凝 
固 结束 后 , 存在 奥 氏 体 相 向 石墨 相 的 转变 . 

表 2 将 温度 冷却 至 740 'C 时 的 模拟 结果 (图 1g) 
中 的 石墨 体积 分 数 和 平均 半径 与 实验 结果 (图 1h) 以 
及 根据 平衡 相 图 的 杠杆 定律 计算 值 进行 了 比较 . 模 
拟 结果 中 的 石墨 分 数 和 平均 半径 采用 程序 计算 和 
Imagetool 软件 2 种 方法 测 得 , 实验 照片 中 的 石墨 体 
积分 数 和 平均 半径 由 Imagetool 软件 测 得 . 从 表 2 可 
以 看 出 , 关于 模拟 的 结果 , 用 程序 计算 和 Imagetool 
则 量 的 数值 非常 接近 . 模拟 的 石墨 体积 分 数 略 低 于 
杠杆 定律 计算 的 结果 . 这 是 因为 杠杆 定律 是 基于 
完全 平衡 的 相 变 过 程 , 而 CA 模拟 包含 了 凝固 过 程 
非 平 衡 的 动力 学 因素 ; 另 一 方面 , 模拟 值 又 略 高 于 
实验 结果 , 这 是 因为 实验 是 一 个 三 维 空间 , 图 1h 区 
域内 的 石墨 不 一 定 为 石墨 球 最 大 的 横 截 面 . 而 模拟 
是 在 二 维 区 域 进行 的 . 综合 这 两 方面 的 比较 , 可 以 
得 出 本 模型 的 模拟 结果 是 合理 的 . 
22 石墨 与 奥 氏 体 相互 作用 机 理 

为 了 研究 凝固 过 程 中 石墨 与 奥 氏 体 之 间 相 互 
作用 与 竞争 生长 的 机 理 , 本 工作 模拟 了 Cu=4.65% 过 
< 唱 球 墨 铸铁 的 凝固 过 程 . 设 定 区 域内 4 个 石墨 和 
一 个 奥 氏 体 枝 晶 的 位 置 , 模拟 区 域 为 200x200 个 网 
格 , 网 格 尺寸 为 2 um. 初始 温度 为 该 成 分 所 对 应 的 
石墨 液 相 线 温 度 1329 "C, 并 以 3 'C/s 的 速率 进行 
冷却 , 直到 温度 降 至 (ZL,-10) C, 随后 设 温度 为 恒 
XE. 图 3 为 模拟 的 石墨 与 奥 氏 体形 貌 的 演化 以 及 深 
度 场 的 变化 . 由 图 3 可 以 看 出 , 在 凝固 开始 阶段 , 当 
温度 低 于 石墨 液 相 线 温度 时 , 石墨 在 液 相 中 形 核 ， 
随 着 石墨 生长 吸收 溶质 C, 其 周围 液 相 的 C 浓 度 低 
于 远离 石墨 的 液 相 区 域 . 并 且 石 墨 I 和 II 周围 的 浓 
度 略 低 于 石墨 1( 图 3a), 3A DS] 778 E TERI TTE E [e] 
消耗 近邻 区 域内 的 溶质 C, 使 得 附近 区 域 的 C 浓 度 


"T 


表 2 模拟 的 石墨 体积 分 数 和 平均 半径 与 杠杆 定律 计算 值 以 及 实验 值 的 比较 


Table 2 Comparisons of simulation, lever rule and experiment for the graphite volume fraction and average radius of a hy- 


pereutectic SG cast iron with Co=4.71% 


Method 


Graphite volume fraction / 96 


Average radius / um 


Simulation result of program calculation 
Simulation result of Imagetool measurement 
Lever rule calculation 


Experimental result of Imagetool measurement" 


14.16 21.24 
14.10 21.19 
14.59 21.55 
11.44 18.52 
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QE AR: 球墨 铸铁 凝固 


BUR. 随 着 温度 降低 至 共 唱 温度 , 奥 氏 体形 核 (图 


3b). 在 奥 氏 体 生 长 过 程 中 先 包围 
墨 , 接着 包围 I 和 II 号 石墨 ， 


与 其 较 近 的 I 号 石 
JE JS II, 并 继续 


生长 . 由 于 奥 氏 体 生长 过 程 中 排出 溶质 C, 而 石墨 的 


生长 吸收 溶质 C, 使 得 


于 石墨 周围 的 浓度 (图 3b~g). 


凝固 过 程 中 石墨 的 生长 速度 和 


变化 如 图 4 所 示 ， 


奥 氏 体 周 


围 区 域 的 C 浓 度 高 


径 随 时 间 的 


时 中 生长 速度 通过 石墨 平均 半径 


的 增长 速度 计算 获得 . 从 图 4a 中 可 以 看 出 ,凝固 开 
台阶 段 , 3 个 石墨 的 生长 速度 都 较 快 . 随 着 石墨 在 液 


相 中 生长 吸收 溶质 C, 使 得 液 相 
降 , 石墨 周围 形成 浓度 梯度 . 根据 式 (7), 石墨 界面 C 


中 的 C 浓 度 快速 下 


浓度 的 降低 导致 石墨 生长 速度 减 慢 , 随后 达到 相对 


稳定 的 生长 速度 ( 
度 基 本 保持 不 变 ( 


图 4a), 这 是 因为 此 时 石墨 周围 浓 


ChinaXiv 合 作 期 刊 


显 微 组 织 的 元 胞 自动 机 模拟 155 


石墨 I 周 于 的 浓度 略 高 于 石墨 I 和 II (图 3a~c), 所 
以 其 生长 速度 略 快 于 其 他 2 个 石墨 . 然后 , 奥 氏 体 开 
始 包 围 石 墨 I( 图 3d), 奥 氏 体 生长 时 排出 的 洲 质 促 
进 了 石墨 1 的 生长 , 使 石墨 I 的 生长 速度 明显 加 快 . 
当 石 墨 1 完 全 被 奥 氏 体 包裹 之 后 , 石墨 1 的 生长 速度 
迅速 减 慢 至 低 于 开 号 石墨 的 生长 速度 , 并 趋 于 一 个 
相对 稳定 值 . 这 是 因为 石墨 在 奥 氏 体 中 生长 时 必须 
通过 奥 氏 体 壳 吸收 洲 质 , C 在 奥 氏 体 中 的 扩散 速度 
较 慢 , 使 得 被 奥 氏 体 完全 包 于 后 的 石墨 1 的 生长 速 
度 减 慢 . 而 此 刻 石墨 开 仍 在 液 相 中 生长 , 所 以 其 生长 
速度 快 于 石墨 I 随 着 凝固 的 进行 , 奥 氏 体 逐 渐 接近 
f 88 ILU TIL, 这 2 个 石墨 在 与 奥 氏 体 接触 和 相互 作 
时 的 生长 速度 变化 规律 与 石墨 I 的 情况 类 似 . 当 
其 被 奥 氏 体 接近 和 包围 时 ， 生长 速度 明显 加 快 被 


—- 


图 3b-d). 从 凝固 开始 到 11 s 之 前 ， 


完全 包围 之 后 , 生长 速度 迅速 减 慢 至 一 个 相对 稳 态 


图 3 模拟 的 C=4.65% 过 共 唱 球墨 铸铁 凝 


司 过 程 中 石墨 与 奥 氏 体 相 互 作用 与 竞争 生长 


Fig.3 Interactive and competitive growth between austenite dendrite and graphite nodules for a hypereutectic SG cast iron 
of Co=4.65% at solidification times of s (a), 6 s (b), 11 s (c), 12 s (d), 13 s (e), 19 s (f), 20 s (g) and 22 s (h) (Num- 


bers in the figures show the local carbon concentration) 
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Fig.4 Time evolution of growth rate (a) and radius (b) of nodules I, II and III in Fig.3 
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值 . 从 图 4b 中 可 以 看 出 , 在 凝固 开始 时 , 3 个 石墨 的 
半径 接近 , 随 着 凝固 的 进行 , 石墨 I 的 半径 大 于 其 他 
2 个 石墨 , 当 3 个 石墨 分 别 被 奥 氏 体 包 围 时 , 其 半径 
变化 规律 与 生长 速率 的 变化 规律 相 类 似 . 此 外 , 在 
整个 凝固 过 程 中 , 石墨 II 的 生长 速度 和 半径 小 于 其 
他 2 个 石墨 , 这 是 因为 该 石墨 受到 周围 2 个 石墨 的 
影响 , 因此 , 最 终 石 墨 II 的 尺寸 最 小 , 并 呈现 出 椭圆 
形 的 形 貌 . 
2.3 冷却 速率 对 球墨 铸铁 凝固 过 程 的 影响 

在 实际 的 凝固 过 程 中 , 对 于 成 分 一 定 的 球墨 铸 
Ek, 不 同 的 冷却 速率 会 影响 石墨 球 的 尺寸 和 尺寸 分 
^g, 从 而 影响 铸件 的 宏观 性 能 . 本 工作 模拟 分 析 了 
冷却 速率 对 Co=4.31% 的 球墨 铸铁 凝固 组 织 中 石墨 
球 尺 寸 的 影响 . 模拟 区 域 和 网 格 尺寸 与 图 3 相同 . 石 
墨 个 数 由 文献 [23] 中 的 实验 测量 数据 确定 , 并 采用 
公式 和 =(NW2.37)” 将 3D 的 石墨 个 数 转换 成 2D 区 
域内 的 个 数 , 其 中 入 和 和 Wy 分 别 为 单位 面积 和 单位 体 
积 的 石墨 个 数 . 此 外 , 文献 [23] 认 为 直径 小 于 5 pm 的 
颗粒 多 为 夹杂 物 和 微 孔 . 因此 , 本 工作 在 从 文献 [23] 的 
数据 计算 石墨 个 数 时 没有 把 直径 小 于 5 hm 范围 内 的 
颗粒 数目 计算 在 内 , 由 此 得 到 的 4 和 1.5 mm 壁 厚 条 件 
下 石墨 球 的 形 核 密度 分 别 约 为 700 和 1400 mm^. 单 
位 面积 奥 氏 体 枝 唱 个 数 NMw 根 据 文 献 [18] 中 形 核 模 型 
公式 N ,, 74.81x10715.33x 10 *(dT/df)-8.7x 10'(dT/dr)" 
计算 得 到 . 由 于 文献 [23] 中 没有 给 出 不 同 壁 厚 所 对 
应 的 冷却 速率 , 本 工作 模拟 过 程 中 的 散热 速率 根据 
文献 [24] 中 的 实验 值 确定 , 4 和 1.5 mm 壁 厚 的 球墨 
铸铁 凝固 时 的 散热 速率 分 别 为 16 和 21 'C/s, 区 域 中 
的 冷却 速率 受到 散热 速率 和 石墨 与 奥 氏 体 析出 释放 
的 潜 热 共同 影响 . 根据 随 着 冷却 速率 增 大 , 石墨 和 奥 
氏 体 形 核 时 的 过 冷 度 增 大 的 实验 现象 "对 于 4 mm 
壁 厚 的 球墨 铸铁 , 模拟 时 初始 温度 设 为 初始 成 分 所 
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对 应 的 石墨 液 相 线 温 度 Tu, 并 伴随 着 初生 石墨 在 
液 相 中 形 核 . 当 区 域 温度 降 至 (us-15) 'C 时 , 奥 氏 体 
形 核 . 当 温 度 进一步 降 至 (Du-35) "CH, 共 唱 石墨 形 
核 . 对 于 1.5 mm 壁 厚 的 球墨 铸铁 , 初生 石墨 、 奥 氏 体 
和 共 唱 石墨 的 形 核 温 度 分 别 为 (Tia-20), (Kxs-25) 和 
(ur-45) °C. 石墨 与 奥 氏 体 的 形 核 位 置 随机 确定 . 图 
5 和 6 分 别 为 Co= 4.31% 的 球墨 铸铁 在 不 同 散热 速率 
条 件 下 石墨 的 尺寸 分 布 和 形 貌 图 . 可 以 看 出 , 随 着 
冷却 速率 增 大 , 石墨 球 个 数 增多 , 石墨 尺寸 分 布 的 
峰值 更 为 明显 , 并 且 峰 值 从 右 向 左 向 着 石墨 半径 较 
小 的 区 域 移动 , 石墨 球 的 最 大 直径 也 有 所 减 小 . 这 
说 明 随 着 冷 速 增加 , 大 多 数 石墨 球 的 半径 减 小 . 当 
冷却 速率 较 小 时 , 石墨 的 尺寸 分 布 呈 现 出 2 个 峰值 ， 
这 是 由 于 初生 石墨 和 共 晶 石墨 的 两 步 形 核 机 理 引 
起 的 . 而 随 着 冷却 速率 的 增 大 , 石墨 两 步 形 核 之 间 
的 时 间 缩 短 , 初生 石墨 在 共 晶 石墨 形成 之 前 的 生长 
TPR, 从 而 形成 单一 峰值 的 尺寸 分 布 特征 . 将 图 5 和 
6 的 模拟 结果 与 文献 [23] 报 道 的 实验 数据 进行 比较 ， 
模拟 和 实验 得 到 的 石墨 尺寸 分 布 与 形 貌 随 着 冷却 
速率 变化 的 规律 相符 合 . 

为 了 进一步 研究 冷却 速率 对 石墨 球 平 均 半 径 的 
影响 , 在 图 5 的 基础 上 增加 了 3 和 2 mm 壁 厚 条 件 下 的 
模拟 . 单位 面积 石墨 与 奥 氏 体 个 数 的 计算 方法 与 图 5 
的 相同 , 3 和 2 mm 壁 厚 时 石墨 形 核 密度 分 别 为 880 和 
1160 mm^, 散热 速率 分 别 取 18 和 20 'C/s. 3 mm RE JE 
条 件 下 , 初生 石墨 、 奥 氏 体 和 共 唱 石墨 的 形 核 温度 
分 别 为 (Tiu-8), (T4,—19):81(7,-39) 'C.2 mm EE JE 
时 , 初生 石墨 、 奥 氏 体 和 共 晶 石墨 的 形 核 温 度 分 别 
为 (Tiu-16), (Tw 一 23) 和 (43) C. 其 他 模拟 条 件 均 
与 图 5 相同 . 图 7 是 模拟 的 石墨 平均 半径 与 凝固 时 
间 的 关系 . 可 以 看 出 , 随 着 凝固 时 间 缩 短 , 石墨 的 平 
均 半 径 减 小 . 因为 凝固 时 间 越 短 , 意味 着 冷却 速率 
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5 模拟 的 Co=4.31% 的 球墨 铸铁 在 740 C 时 石墨 尺寸 分 布 随 散热 速率 的 变化 
Fig.5 Simulated nodule size distribution for a SG iron of Co=4.31% at 740 °C with extraction rates of 16 'C/s for 4 mm 


plate (a) and 21 'C/s for 1.5 mm plate (b) 
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图 6 模拟 的 C=4.31% 的 球墨 铸铁 在 740 'C 时 显 微 组 织 形 貌 随 散热 速率 的 变化 


Fig.6 Simulated microstructures for a SG cast iron of Co=4.31% with heat extraction rates of 16 °C/s (a) and 21 'C/s (b) 
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球墨 铸铁 石墨 平均 半径 与 凝固 时 间 的 关系 


Fig.7 Comparison of the average graphite radius as a func- 


tion of solidification time between simulation (C= 
4.3196) and experiment (Fe-3.57%C-2.64%Si)™ re- 
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li 


越 大 , 形成 的 石墨 个 数 越 多 , 而 对 于 初始 成 分 一 定 
的 球墨 铸铁 , 凝固 结束 时 石墨 的 体积 分 数 变 化 不 
X. 因此 , 冷却 速率 越 大 , 最 终 的 石墨 球 平均 半径 越 
小 . 从 图 7 还 可 以 看 出 , 模拟 结果 与 实验 结果 所 显示 
的 规律 相 一 致 . 


3 结论 
(1) 改进 了 前 期 工作 建立 的 多 相 元 胞 自动 机 


(multi-phase cellular automaton, MCA) 模 型 , 模拟 以 
离异 共 唱 方式 凝固 的 球墨 铸铁 的 显 微 组 织 演 化 . 在 
模型 中 采用 基于 局 部 溶质 平衡 的 方法 计算 石墨 与 
奥 氏 体 的 生长 , 将 界面 平衡 成 分 与 实际 成 分 的 差 值 
作为 生长 的 驱动 力 . 并 且 在 计算 石墨 的 生长 动力 学 
中 考虑 了 石墨 与 Fe 的 密度 比 . 应 用 该 模型 可 以 模拟 
出 与 实验 观察 吻合 良好 的 离异 共 章 球墨 铸铁 凝固 


显 微 组 织 形 貌 . 模拟 的 石墨 体积 分 数 和 平均 半径 比 
实验 测量 结果 略 高 ` 而 比 杠杆 定律 计算 值 略 低 , 体 
现 了 模拟 结果 的 动力 学 效应 . 

(2) 模拟 分 析 了 球墨 铸铁 凝固 过 程 中 石墨 与 奥 
氏 体 两 相 之 间 相 互 作 用 和 竞争 生长 的 机 理 . 结果 表 
8j, 奥 氏 体 在 液 相 中 的 生长 可 以 促进 附近 石墨 的 后 
T, 两 相 的 生长 受 C 扩 散 控 制 , 当 石墨 被 奥 氏 体 包 围 
之 后 , 由 于 C 在 奥 氏 体 中 的 扩散 系数 远 小 于 其 在 液 
相 中 的 扩散 系数 , 石墨 生长 速度 明显 减 慢 . 相 邻 的 
墨 之 间 也 存在 相互 竞争 生长 的 现象 . 此 外 , 奥 氏 
体 与 石墨 球 接触 以 前 以 校 晶 的 形态 生长 ; 当 奥 氏 体 
包围 石墨 球 之 后 , 以 团 状 方式 生长 而 失去 原先 的 校 
蝇 形 态 . 

(3) 模拟 分 析 了 冷却 速率 对 凝固 结束 时 石墨 球 
的 尺寸 分 布 以 及 平均 半径 的 影响 . 结果 表明 , 当 冷 
却 速率 较 小 时 , 石墨 尺寸 分 布 呈现 2 个 峰值 , 这 是 由 
于 初生 石墨 和 共 品 石墨 两 步 形 核 所 产生 . 随 着 冷却 
速率 的 增 大 , 石墨 两 步 形 核 之 间 的 时 间 缩 短 , 初生 
墨 在 共 唱 石墨 形成 之 前 的 生长 有 限 , 石墨 尺寸 分 
布 变 为 一 个 峰值 , 并 且 峰 值 向 着 半径 较 小 的 区 域 移 
动 . 冷却 速率 越 大 , 凝固 时 间 越 短 , 石墨 的 平均 半径 
越 小 . 模拟 结果 与 文献 中 的 实验 结果 吻合 良好 . 

感谢 美国 Ohio State University 的 Doru M. Stefanescu 教授 关于 本 
工作 的 有 益 讨论 . 
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